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£\j , Abstract 

Compaction in reactive porous media is modelled as a reaction-diffusion process with a mov- 
ing boundary. Asymptotic analysis is used to find solutions for the coupled nonlinear compaction 
equations, and a traveling wave solution is obtained above the reaction zone. 
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^ ■ 1 Introduction 

The accurate modelling of compaction and reactive flow in porous media such as sand and shale 
sediments is very important in civil engineering and oil industry [1] . The general mathematical model 
of compaction and reactive mineral flow considers the fluid-sediment system as a porous medium 
consisting of multiple mineral species. The interstitial volume of the porous solid phase is saturated 
with pore fluids. Due to the action of gravitational overburden loading and the density difference 
between the two phases, the solid phase compacts by reducing the porosity (volume fraction of the 
pores) , thus leading to the expulsion of the pore fluid out of porous media. During compaction and 
continuous burial, the mineral species such as water-rich clay smectite react with pore fluids and 
£T) ■ are then transformed into a more stable mineral species such as illite, releasing free water into the 

porous environment [1,2]. In this paper, a reaction-diffusion model together with some asymptotic 
analysis is presented. 



2 Mathematical Model 

Let the volume fractions of solid reactant species (smectite) and water be ip, 0, respectively. By 
proper non-dimensionalization and appropriate scalings [3,4,5] in a 1-D basin < z < h(t), where 
h(t) is the ocean floor and z = is the basement rock, we can write down the non-dimensional 
compaction equations as 

It = -eP< h -*-*">ii> - r\m^) m (^ - <P)h- (i) 

1-00 00 

4>t =*[(-£■ r(0 z -0)] z + ^V('™*V, m>7. (2) 

00 P 



The boundary conditions are 



Z - = 0, at 2 = 0, (3) 

0=00, ,1p=1p , 

h(t) = s + ^^(±r^ z - 0) at z = h(t), (4) 

1 -00 00 
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where A = O(l) and > 1 are compaction constant and non-dimensional activation energy of the 
one-step dehydration mineral reaction with a critical temperature z* and the amount a — 0(1) of 
free water released from the reaction, s is the rate of new sediment accumulation at the basin top, 
and thus can be taken as a prescribed constant (s = 1). <fio and ipo are the initial values of and tp at 
the ocean floor z = h(t), respectively. In addition, all the variables {cf), tjj, z, t) and the parameters 
(A, m, (fto, ip , s) are non-negative. The reaction term exp[/3(/i — z — z*)] is only dominant within a 
thin region of a width of 0(1/ /?) near h — z — z* « 0, in other words, the reaction region is located 
at z « 9* defined as 

9* = h-z*. (5) 

Clearly, the present problem is a non-linear diffusion problem with a boundary h(t) moving at a speed 
of h(t), which can be solved numerically by using the predictor /corrector implicit finite-difference 
method. 



3 Travelling Wave Solution 

Numerical simulations [3,6] and real data observations [7] imply that the moving boundary h(t) 
moving essentially at a nearly constant speed h = c, although the specific value c depends on the 
boundary conditions, and is yet to be determined. The constant moving boundary implies there 
exist travelling wave solutions. Thus, we define a new variable by 

( = z-9* = z-h(t) + z*, -6*<(<z*, (6) 

so that the model equations (1) and (2) become 

-c^c = -e _/3C V> - T^-[^(^) m (^ ~ 0)k- (7) 

f - 00 00 

- c ^ = A[(-fr(^-0)] f + ^ e -^v, (8) 

0o P 

The fact that /3 ~ m ^> 1 allows us to seek asymptotic solutions in different regions with these 
distinguished limits. /? > 1 implies that the reaction is essentially restraint in a very narrow zone 
with a width of 0(l//3) at z — 9*. Above this region (z > 9*), we have £ > so that cxp(-^C) -> 0. 
Below this reaction region (( < 0), the reaction is essentially completed (i.e., the volume fraction of 
smectite — > 0) and consequently ip exp(—f3() — > 0. We shall see that this is true below in equation 
(24) because ip — > as ( — > — oo (or r) — > — oo). 



Outer Solutions 

In the outer region above the reaction zone (£ > 0), the reaction terms are negligible, then we have 

cA = T ^rM^-) m (<t><-<f>)h (9) 
J- - <Po <Po 

-c0 c = A[(|-r(^-0)] c . (10) 
The integrations together with top boundary conditions (4) give 

and 

c<P + A(-f- ) m (0 c - 0) = c0o + (c - *)(1 - 0o), (12) 

whose solution can be further written as a quadrature. This solution is only valid in the region 
above the reaction zone (z = 6*). On the other hand, the travelling solution will not be appropriate 
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in the region below the reaction zone (£ < 0) because the exponential term (4>/(j) ) m <C 1 due to 
4> < (f>o and m 3> 1. However, we can define a typical value of (f>* by 

<j>* = foe^ , (13) 

so that <f> ~ 0* in the region below the reaction zone (£ < 0). The reaction term is still negligible. 
Thus we rewrite equations (1) and (2) in terms of $ defined by 

— _ — ln _ 

^ — ^ em — ^ oe m , (14) 

then we have 

* - ~^M^i*' - ')'■■ <15 » 

ft, = V[ t '(i»,-l)!„ (16) 
By using 1/m <C 1, the above equations becomes 

& » 0, (17) 

$ t + Ae*$ z = 0. (18) 

It is straightforward to write down the solution together with the boundary condition $ z = in at 
z = 0. We have 

, . 1 + m,z , ,„ „. 

^IT^a^' (19) 

As z = h(t), we have $oo = ln{(l + mh)/(l + mXt)} — > ln(h/Xt) = ln(c/A) as t — > oo or mh(t) — > oc 
due to dh/dt = c = const and h = at t = 0. However, z = /i(£) is not usually reached since solution 
(22) is below the reaction region. When z 1, the above solution shall match the inner solutions 
as n — > — oo. 

Inner Solutions 

In the reaction zone, we use the stretched variables defined by 

77 = /3C + ln/3, = 0*e^ =</, oe ^r la , (20) 
so that equations (7) and (8) become 

-c^ = -e-"V - ^^y[tfe*(A*, - 1)]„. (21) 

-c4>*% = \<p[e*(A*„ - 1)]„ + ^e- r V, (22) 
where ^4 = Jj = 0(1). By using <C 1, equation (21) becomes 

= e-"V, (23) 

whose solution is 

V> = Cexp[-i e - 7 '], C = Vocxp[ie^ z *- ln ' 3 )], (24) 

where C depends on c, and c will be determined later in (29) by matching. Substituting this solution 
into (22), integrating once from — oo to -q and using $ — > $oo as r? ^ — oo, we get 

c<j>*$ + \<f>*e*{A<S> TI -l)-B = -^C, (25) 



where B = c</>*$ 00 — \<f>* exp(<I> 00 ). As 77 — » 00, we have 



[of* + \<fe*{A* n - l)]t~ = --^C, (26) 



which implies a jump through the reaction region. 
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Matching 

By rewriting solution (25) in terms of <f), we have approximately 

c<f> + X(^) m ^c - 4>) » c0*$oo - A0*e*~ - C^, (27) 
00 A 

and matching this to solution (12), we have 

af>o + (c - s)(l - 0o) = c4*<S>oo - A0*e*°° - C^, (28) 
which determines c, leading to 

»Mi) 

l-^oo + ^Z (l-0*$oo + ^)' 

Since $oo is a function of c, we now have an implicit equation for c, which depends essentially on 
the initial values of 0o and Vo- 

In summary, although it is very difficult to seek directly solutions for the couple nonlinear 
reaction-diffusion equations, we use a hybrid method to get the asymptotic solutions and travel- 
ling wave solutions in different regions. The matching of these solutions can thus determine the 
boundary moving velocity h(t) — c, which shows how the reaction inside the porous media affect 
the evolution of its top boundary. 
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